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We introduce a primary thermometer which measures the temperature of a Bose-Einstein Condensate in 
the sub-nK regime. We show, using quantum Fisher information, that the precision of our technique 
improves the state-of-the-art in thermometry in the sub-nK regime. The temperature of the condensate is 
mapped onto the quantum phase of an atomic dot that interacts with the system for short times. We show 
that the highest precision is achieved when the phase is dynamical rather than geometric and when it is 
detected through Ramsey interferometry. Standard techniques to determine the temperature of a 
condensate involve an indirect estimation through mean particle velocities made after releasing the 
condensate. In contrast to these destructive measurements, our method involves a negligible disturbance of 
the system. 

Temperature is a crucial concept in quantum physics. Paradigmatic phenomena such as superconductivity 
and Bose-Einstein condensation only occur below a critical temperature. Bose Einstein Condensates 
(BECs) 1 allow the study of quantum effects in systems consisting of up to 10 s atoms by cooling them to 
regimes in which the individual atomic wavefunctions overlap. In this case, the system exhibits quantum beha- 
viour at mesoscopic scales. 

Very low temperatures are also required to observe quantum field theory effects e.g. Unruh-Hawking radi- 
ation 2 . Quantum field theory predicts that the dynamics of spacetime or the presence of horizons produce 
quantum particles from vacuum fluctuations. For example, in the dynamical Casimir effect, a moving boundary 
condition gives rise to vacuum excitations. This effect was recently demonstrated in superconducting circuits 3 . 
Currently several experimental groups in the field of analogue gravity attempt to demonstrate cosmological 
particle creation, dynamical Casimir effect and Hawking radiation in analogue spacetimes produced in BECs 4 . 
However, in a BEC a thermal background is always present due to unavoidable atomic collisions. In these 
experiments, it is crucial to work at low temperatures so that the quantum particles created by effects of emergent 
spacetimes can be distinguished from this thermal noise. Recent experiments to demonstrate the dynamical 
Casimir effect in a BEC 4 , have not been able to ensure that the excitations produced are indeed quantum even at 
temperatures as low as 200 nK. A new method to accurately determine temperatures in the nK and sub-nK 
regime would be of great benefit to these experiments. 

The temperature of a BEC is commonly estimated indirectly by comparing the density profile of the atoms with 
a velocity distribution 1,5-7 . This density profile is determined by absorption imaging which is a destructive method 
that involves releasing the condensate. Although accuracies of ~ 1 % can be achieved in the 100 nK regime 8 , the 
relative error in the measurement of temperatures considerably grows at very low temperatures. The best 
accuracies that have been reported in the nK and sub-nK regime are of ~20% 9 . Several alternative methods to 
measure the temperature of a BEC have been proposed. For instance, phase-contrast imaging 10 can be used to 
determine the temperature through the phase shift of a probe laser beam that interacts with the sample. However, 
at temperatures of a few nK and below this method has low spatial resolution and can only be applied reliably if the 
condensate is allowed to freely expand for a given time 11 . Another example is noise thermometry with two 
coupled Bose-Einstein condensates 12 , which has been implemented in a regime in which quantum fluctuations 
are small (50-80 nK). Therefore, absorption imaging is the standard technique for measuring temperatures in the 
nK regime 8 and, in general, for probing the condensate 13 . 

Recent theoretical and experimental developments in cold atom gases have lead to the experimental demon- 
stration of systems consisting of mixtures of two different atomic species or two hyperfine states of the same 
species 1418 . In particular, it is in principle possible to create an atomic quantum dot by trapping a few atoms in a 
tight trap provided by a laser beam 20 26 . The number of particles in the dot is determined by the interatomic 
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interactions which can be modulated experimentally, through a 
Feshbach resonance 1719 . In this paper we show how an atomic 
quantum dot immersed within the condensate can be used to mea- 
sure the temperature of the condensate with high precision and a 
negligible disturbance of the system. The scheme can be implemen- 
ted using an optical lattice with a single atom per site 27 , where each 
atom interacts at a different time with the condensate. This system 
constitutes a primary thermometer which measures the tempera- 
tures ranging from the sub nK regime to the condensate critical 
temperature. 

The outline of the paper is as follows. We start by showing how a 
two-level system interacting with a quantum field in a thermal state 
through the standard Jaynes-Cummings Hamiltonian can acquire a 
dynamical phase that depends on the temperature of the field. This 
phase is then read out through Ramsey interferometry. We use the 
Quantum Fisher Information to quantify the accuracy in the tem- 
perature estimation. We show that with realistic experimental para- 
meters, the precision of our thermometer improves on state of the art 
thermometry in the nK and sub nK regimes already after a small 
number of measurements. Since each dot embedded within the con- 
densate provides one measurement, it is important to ensure the 
statistical independence of the individual measurements, by making 
the dots interact independently with the condensate. By preparing 
the dots in an entangled state, it is possible to further improve the 
precision and reach the Heisenberg limit 28 , that is the limit imposed 
fundamentally by quantum theory. Finally, we present details of 
possible physical implementations of the thermometer using atomic 
quantum dots coupled to phonons in a BEC (see Fig. 1). 

While the idea of probing several aspects of the condensate with 
different two-level systems and interactions has been introduced in 
earlier works (see, for instance 29,30 ,) our work represents a step forward 
in several directions, which we outline here. We show that with real- 
istic experimental parameters the system can be described by a Jaynes- 
Cummings model, where the dots are effectively coupled to a single 
mode of the quantum phase fluctuations of the condensate. Moreover, 
we show that is possible to achieve a regime in which spontaneous 
transitions of the dots are negligible (see Methods section) and thus 



the time evolution of the system is given simply by dynamical 
and geometrical phases. We study both phases and analyse the 
advantages of using Ramsey vs. Mach-Zehnder interferometry to 
estimate relative phases. We use Quantum Metrology techniques 
to quantify the precision in each scenario. We show that the best 
precision is achieved when using dynamical phases and Ramsey 
interferometry. Furthermore, we show that this precision can 
improve on standard methods to determine temperature in the nK 
and sub nK regimes. 

Results 

We now proceed to introduce our methods and results. We consider 
the Jaynes-Cummings model that describes the interaction of a two- 
level system of frequency gap Q. d and one mode of a quantum field 
with frequency Q. a 

H; C {t)=hQ a a ! a + h^a z + hg{a-a < e i0 + a + ae- i0 ), (1) 

where a, a 1 are creation and annihilation field operators, h = h/(2n) 
is the reduced Planck constant, a~ = (1/2)((T X ± ia y ) and a x and a y 
are Pauli matrices. The coupling strength between the field and atom 
is given by g, the detuning is a = |Q„ — Q. d \ and the phase 0 = kx— 5t 
is a function of time and the position x of the atom. As we will explain 
in more detail later, with suitable boundary conditions an atomic dot 
interacting with the phonon field of a BEC can be described with this 
simple interaction in the nK and sub-nK regime. 

The temperature of the field will be determined through the phases 
acquired by the atom through its evolution under the Hamiltonian 
(1), assuming the adiabatic approximation (see Methods). The total 
accumulated phase y can be split into a geometric and a dynamical 
part. It has been shown that the geometric phase can be used to gain 
information about the quantum state of a bosonic scalar field. For 
pure field states, the phase encodes information about the number of 
particles in the field 31 . In particular, for initially squeezed states, the 
phase depends also on the squeezing strength 32 . More recently, the 
geometric phase was employed to estimate the temperature of a field 
inside a cavity 33 . We will apply Quantum Fisher Information to show 
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Figure 1 | (a) Sketch of the experimental setup: several atomic quantum dots are embedded within a BEC reservoir. They are coupled at different 
times through a Raman transition to the phononic fluctuations of the BEC. The use of many dots allows one to implement many measurements and 
improve the accuracy, (b) Ramsey interferometry scheme to measure the relative dynamical phase and hence measure temperature, (c) Ramsey 
interferometry scheme with entangled input states. Entanglement can in priciple improve the precision, reaching the Heisenberg limit, (d) Mach-Zehnder 
interferometer scheme with entangled input states. (A). Input state, (B). Interaction switched on, (C). Final state (D). Readout. 
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here that estimating the dynamical phase is the optimal way for 
measuring temperatures. We consider the case where the field is 
initially in a thermal state p F = p„|n)(w| with temperature T. 
The probability distribution is given by p n = (e~ F )" (1 — e _F ) where 

F = - — — , and kn is Boltzmann's constant. Here we have used the 
k B T 

number state basis spanned by | «) = (a 1 )" |0), where n is the number 
of phonons of frequency fl a . We find (see Methods), that the geo- 
metric phase is 



5 l {e F -l) 



and the dynamical phase 
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As explained in Methods, in order to compare the phases we consider 
t = 2n/S, since cyclic evolution is necessary to compute the geometric 
phase given in Eq. (2). We find that the dynamical phase is always 
larger than its geometric counterpart 

y_c_ 



(4) 



In Fig. 2a) we plot the dynamical and geometric phases for realistic 
experimental parameters. One observes that for temperatures in the 
sub nK regime, the geometric phase is negligible in comparison to the 
dynamical phase. Furthermore, we will show that the precision 
obtained through the dynamical phase in the measurement of the 
temperature is higher than the precision attainable with the geomet- 
ric one. 

Since the phases acquired by the system are global and cannot be 
measured directly, it is necessary to setup an interferometric scheme 
that will allow us to measure relative phases. In particular, the geo- 
metric phase can only be detected using a Mach-Zehnder interfero- 
meter. In this type of interferometer the qubit is split into two 
different trajectories and is made to interact with a thermal field in 
one arm of the interferometer acquiring a relative phase that can be 
measured at the interferometer's output. Now we will use standard 
techniques of phase estimation in Quantum Metrology to estimate 
the precision of the temperature measurement in this experiment 
(see Methods) and determine which type of phase (geometrical or 
dynamical) is more convenient. 

We consider that in one arm of the Mach-Zehnder interferometer 
the qubit, initially in state \g), interacts with a field in state p T . In the 
second arm we assume that the qubit in state |e) does not interact 



with the field. Since, up to the first order in g/8 the state 
\s) (s\ ® Pt — ^2 Pn\n—){n—\, the relative phase at the output will 
therefore be given by y~ = y^ + yT. . For the sake of comparison, we 
first assume that the dynamical phase can be cancelled 34 and we 
calculate the quantum Fisher information to estimate the precision 
on the temperature measurement through the geometrical phase. 
Since in this case F(p T ) = {pry^ ) 2 , the Cramer-Rao bound in the 
geometric case employing Mach-Zehnder interferometry, denoted 
by <5T°, yields 

i r<5V-2) 



Md T y G 



MFg 2 



(5) 



Here the bound is obtained using a series expansion in g/S. Following 
the same procedure we find bounds to the quantum Fisher informa- 
tion for the case where the temperature is measured using the 
dynamical phase. We find that: 
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From Eq. (5) and (6) it follows that the error in estimating the 
temperature using the dynamical phase is smaller ST^kST^. 
Furthermore, using the dynamical phase allows us to simplify the 
setup by using Ramsey rather than Mach-Zehnder interferometry. 
Ramsey interferometry does not involve spatial splitting the atomic 
wave function to undergo different trajectories and therefore is less 
demanding from the experimental viewpoint. 

In Ramsey interferometry the two-level system is prepared in a 
superposition state, |f) = \g) + je), and is allowed to interact 
with the quantum field in a thermal state px during a time t. 
After the interaction, the probability of finding the two-level 
system in the excited state depends on the dynamical relative 
phase, r D , picked up by it. Under weak adiabatic evolution the state 

IT)(TI®Pr- X/"(l"~) + l" + ))(( n ~l + < n + l) aa l uires the 
relative phase, To = y£ — y^ , which is very sensitive to the temper- 
ature in an ultralow regime such as nK and sub-nK. Here, y£ , y^ are 
the dynamical phases acquired by |+) and |— ) respectively (see 
Methods section). The Fisher information in this case is given by 

F(Pt) = (StYd ~ ^tYd ) 2 = 4 (drYo Y ■ Therefore, the achievable 
precision in the temperature measurement is increased by a factor 
of 2 for Ramsey interferometry, i.e. ST^ = (1/2) ST^. 

In Fig. 2b we plot the dependence of T D on the detuning <5 and the 
coupling strength g for the relevant experimental regime that we 
consider below. Note that the phase is sensitive to gib and is essen- 
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Figure 2 | (a) Comparison of geometrical (dahed, blue) vs. dynamical y^ (solid, red) with Sl a = 2 n X 10 Hz, g = 2 n X 0.2 Hz, d = 2 n X 2 Hz and 
c = 5 mm/s. ~8k y^ . The dynamical phase is much more sensitive to the temperature in the sub- nK regime, (b) Dependence of the dynamical 
phase T D on g and <5 at T= 0.5 nK. T D is sensitive to the ratio gl 6. (c) T D (green, solid) and relative error tST/Tafter 1000 (blue, dashed), 3000 (red, dotted) 
and 10000 (yellow, dash-dotted) measurements. The parameters are the same as in (a). The precision in the best case is around 1% and can be improved by 
increasing the number of measurements and using entangled input states. 
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tially determined by this ratio and Sl a . In Fig. 2c we plot the relative 
error ST^/T after 1000, 3000 and 5000 measurements. In the best 
scenario, the relative error achieved is around 1% in the range 0.3- 
1 nK, which improves the 20% of the current state of the art ther- 
mometry in the sub nK regime 9 . 

Experimental Implementation 

We now consider a realistic implementation of the above results. We 
will consider a system consisting of a BEC superfluid reservoir in a 
shallow confining trap interacting with an array of atomic quantum 
dots 20,35 . Each dot is created by applying a localised steep potential 
that traps atoms of a different hyperfine state of the same atomic 
species as the BEC. As we will see below this system obeys the Jaynes- 
Cummings Hamiltonian where the coupling is mediated by a Raman 
transition. We would like to point out that two-level systems can also 
be implemented using a deep double-well optical lattice loaded with 
cold atoms in a single-occupancy regime 36 . This system also obeys 
the JC Hamiltonian. However, in this paper we will discuss in more 
detail the atomic dot implementation 20,35 . 

In the absence of atomic collisions, the BEC can in principle reach 
absolute zero temperature. Under this assumption, the condensate is 
commonly described by a classical density function. In realistic scen- 
arios, collisions are always present and therefore, in the superfluid 
regime, the condensate is better described by a mean field classical 
background plus quantum fluctuations. The fluctuations, for fre- 
quencies lower than the so-called healing length, are given by a 
phononic quantum field operator n(x), which can be expanded in 
terms of Bogoliubov modes. Now we proceed to describe the 
quantum dots and their interaction with the phononic field. Each 
dot is created by applying a localised steep potential which traps 
atoms of a different hyperfine state of the same atomic species as 
the BEC. Alternatively, in 36 the two-level systems are created by 
loading cold atoms of a different species in an optical lattice consist- 
ing of double well potentials where only single occupancy of each site 
is allowed. Focusing on the former, we choose a large collisional 
interaction strength g aa inside each dot well, to guarantee that the 
occupation number in each site is either 0 or 1, giving rise to an array 
of two level systems. We also assume that the wells are deep and 
separated enough to neglect direct interaction between different 
sites 35 . 

The atoms in the BEC and the dots can be coupled through a 
Raman transition using external lasers, giving rise to a phonon- 
mediated interaction 20 . The effective Hamiltonian obtained under 
the above assumptions can be written as 20-35 : H = H A + H B + 
H AB . Here H A is the free Hamiltonian of an array of M dots, which 

is given by H A = . _ 1 hCld<j' z /2. The frequency is a function of 
the effective Rabi frequency of the Raman transition Q and the 
number of atoms in the condensate 20 . H B is the free Hamiltonian 
for the phonons H B = hco^a^a^, where io k = c|k| and c is the 
speed of sound. Finally, the interaction Hamiltonian is given by: 

M M 

H AB =-U'Y J < + Kgab -gbb) £ 11 ( X 0 < ( 7 ) 
1 = 1 i=l 

where x' is the position of the dot in the condensate. The coupling 
strengths between dot-condensate and condensate-condensate 
atoms are given by g ab and g bb , respectively. The frequency 5' is a 
function of g a y g aa , the effective Rabi frequency and the detuning. We 
assume that all the coupling strengths g ab , g aa and g hh can be tuned 
with Feshbach resonances. In particular, we will consider that the 
strengths are tuned such that 8' = 0. In this case the system is 
described by the spin-boson Hamiltonian 20 , in which the structure 
of the reservoir is characterised by the spectral density. The ubiquity 
of physical systems that can be described by the spin-boson 



Hamiltonian directly illustrates that the thermometry technique we 
describe can be applied to a range of different physical systems. 

With suitable boundary conditions for the condensate trap -e. g. 
hard-wall or box potentials 37,38 -, the energy separation of the pho- 
nonic modes can be large enough to ensure that each two-level sys- 
tem is only effectively coupled to the mode with closest frequency to 
£lj. We denote this frequency by Q a and apply the Rotating wave 
Approximation (RWA). Notice however that we consider the JC 
model only for the sake of mathematical simplicity. Indeed, if the 
interaction is given by the more general spin-boson model (Eq. (7)) 
our results can be extended by using perturbation theory techniques 
with the multimode field. Assuming that each dot interacts with the 
condensate at different times, the Hamiltonian for each dot is given 
by Eq. (1), where now the dot-field coupling strength g is given by: 



ck 



2hVg bl 



gab -gbb), 



(8) 



where V is the condensate volume. The RWA holds for weak cou- 
plings and small detunings {g,S} « e where t = Q. a + Q^. Taking 
typical values for the condensate length L = 500 [an and the speed 
of sound c = 5 mm/5 and assuming h ^ L, the frequency of the 
phonons is given by Si a = 2n X 10 Hz. With experimental data 
for the Feshbach resonances 1719 the values of g can be tuned within 
a broad range from 2ti X 0.1 Hz to 2n X 10 Hz. Taking g around 
0.1 Hz and a detuning of 2n X 1-2 Hz we are consistent with both 
RWA and the adiabatic approximation. This is the range of para- 
meters explored in the plots. 

Dissipation in this model is characterized by 20 a = —l\ = 



V mc 2 

Y 

1 ] . As can be seen in Eq. 8, the weak coupling regime that 

we are considering entails that g ab is slightly larger than g bb . In par- 
ticular, for g = 0.1 Hz, gablgbb = 1 + 1/70 and then a = 2 • 10~ 4 . 
Therefore, we are in the regime a « 1, where the damping rate F of 
the Rabi oscillations is negligible as compared with the frequency of 
the dot Qj, and the dynamics is well described by perfect Rabi oscilla- 
tions, as in the Jaynes-Cummings model. 

The description of the experiment is the following. The atomic 
dots are prepared in a separable superposition of their internal levels, 
®i=i...m I tf)- They are then coupled to the phonons through Raman 
transitions and allowed to interact with the BEC for some time, 
which in our plots is of the order of 0.1 s. The times can be made 
significantly shorter by relaxing the cyclicity condition, which we 
imposed in order to compare the dynamical phase analysis with 
the geometrical phase analysis. However, it is no longer necessary 
in this setting. The probabilities of excitation of all the dots are then 
measured using standard imaging techniques 27 . In order to ensure 
that the measurements are completely independent, one possibility is 
that each dot interact at different times with the condensate. 
Moreover, the separation in time between interactions should be 
large enough to rule out temporal correlations between the occu- 
pation numbers of the considered mode. These temporal correlations 
drop to 0 after a few ms for typical parameters 39 . Note that this 
implies that the disturbance of the condensate in each measurement 
is negligible. The sequential interactions can be achieved using the 
dependence of the coupling strength on the value of the magnetic 
field near a Feshbach resonance. For the dot that is interacting with 
the condensate, the value of g ab is slightly different from g bb , giving 
rise to a finite value of g. For the rest of the dots, g ab = g bb and g = 0. 
Wheng = 6 = 0, as discussed in 20 , an extra term in the Hamiltonian - 
which is negligible when the interaction is on- becomes relevant. This 
term does not generate transitions in the qubit basis, since it is just a 
dephasing interaction. The ratio between the strength of this dephas- 



ing and the dot frequency is 



hk 



8p b cVm 



:10~ 



. meaning that deco- 
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herence of the dots can be neglected and they can safely be assumed 
to remain in the initial superposition state while the other dots are 
being sequentially let to interact and measured. 

Another way to ensure that different measurements are not cor- 
related is to couple each dot to a different mode of the field. In this 
case, the interactions can be simultaneous, since each dot would 
probe a different thermal distribution with the same temperature. 
Alternatively, an array of dots with different gaps could be used, for 
instance taking advantage of a non-homogenous magnetic field. 
Other sources of imperfections- e.g. laser noise- would also change 
the details of the dot-condensate interaction -in this case, the coup- 
ling strength-, ensuring the statistical independence of the measure- 
ments. The only limitation to the number of probes is that the spatial 
separation must be larger than the healing length of the condensate in 
order to ensure no direct interaction between them 35 . Moreover, 
indirect interaction due to phonon exchange is proportional to g 2 / 
8, and therefore is negligible for times t = 1 18 and g/8 <C 1 . Taking for 
instance a realistic separation between the dots of 300 nm and a 
condensate length of 500 ;im, we can load 1500 impurities in a 1- 
D optical lattice. This number would be further increased in2-D or 3- 
D setups. Since the lifetime of the condensate can be larger than 100 s 
it seems feasible to perform at least a few measurements on each dot. 
Therefore, as can be seen in Fig. 2c) the relative error is well below the 
20% state of the art thermometry in the sub nK regime 9 , and ultra- 
high precisions below 1% can be achieved. Moreover, as explained 
above, by preparing several dots in an entangled input state 22,40 the 
precision would be significantly improved reaching the Heisenberg 
limit. The initial entangled state can be prepared by using a value of 
the magnetic field that strongly couples the dots to the field 41 . Then 
the magnetic field can be varied to decouple the dots and start the 
experiment, as explained above. Notice that, once the dots are 
decoupled, the local dephasing interaction cannot change the non- 
local correlations i. e. entanglement. 

Discussion 

Our minimally invasive technique is suitable to measure with high 
accuracy the very low temperature effects generated in analogue 
gravity scenarios. The Unruh temperature can in principle be mea- 
sured in our setting by accelerating the dots. The idea of demonstrat- 
ing the Unruh effect in BECs was suggested with a different setup 
in 42 . Moreover, the use of a geometric quantum phase as an estimator 
of the Unruh temperature was proposed using a challenging experi- 
mental setup involving atom Mach-Zehnder interferometry 43 . We 
are currently working on implementing our technique to measure 
the Unruh effect in a BEC with higher precision. If a is the accelera- 
tion, the Unruh temperature is given by Tu = h a/(2n c k B ) which in 
our setup give rise to Tu = 2.4 nK for an acceleration ofg -where g is 
the acceleration of gravity on the surface of the Earth. Accelerations 
of several g were achieved in recent experiments 27 . Therefore, with 
our method the detection of the Unruh temperature of the phononic 
BEC bath is clearly within reach of current technology. 

We have shown that the quantum dynamical phase acquired by a 
two-level system interacting with a quantum field in a thermal state 
under the Jaynes-Cummnigs Hamiltonian can be used to measure 
the temperature of the field. Since the Jaynes-Cummings model suc- 
cessfully describes a number of different quantum systems, our ideas 
can be applied to measure the temperature in a great variety of 
physical setups, such as cavity QED, ion traps or circuit QED in a 
weak coupling regime. As a specially interesting implementation we 
choose to present in this paper a cold atom setup. In particular, an 
atomic quantum dot is coupled to the phononic quantum fluctua- 
tions of a BEC. We compute the dynamical phase acquired by several 
independent dots in a Ramsey interferometry scheme and estimate 
the precision by means of the Quantum Fisher Information. We 
show that the phase is sensitive to temperatures in the sub nK regime 
and that the precision improves on the state of the art thermometry at 



those temperatures after a number of measurements that can be 
achieved with current technology. Moreover, the precision can in 
principle be further improved by preparing entangled input states. 

The proposal to infer the temperature of the condensate from the 
phase acquired by an interacting probe, is not restricted to the par- 
ticular BEC model we present here. We chose as a particular example 
a probe coupled to a single mode of a weekly interacting BEC because 
of its mathematical simplicity and elegance. In principle, our ideas 
can also be applied to different cold atom systems, such as Fermi 
gases -where they would be useful in the analysis of magnetic phase 
transitions. Moreover, other states of the field might be considered 
and therefore the readout of the qubit's phase could be a tool to 
analyse the dynamics of a BEC and other cold-atom systems. 

In summary, we have introduced a non-invasive, accurate and 
quantum way of measuring ultralow temperatures in a Bose- 
Einstein Condensate. In addition to improving on state-of the art 
thermometry in the sub nK regime, these techniques provide an 
experimentally feasible method of measuring very-low temperatures 
arising in analogue gravity scenarios, paving the way for the detection 
of analogue Unruh and Hawking temperatures in BECs. 

Methods 

Dynamical and geometrical phases. The time evolution of the system can be found 
by writing H, c (t) = U(B)H^ c U(Bf where U(B) = e i0 "' a andH° c = H ;c (x = 0, t = 0). 
The eigenstates of H® c , are given by n — ) — cos ^— J (n + l)g) — sin ^— J | ne) and 

\n + ) — sin^-^j + + cos (^) \ ne ) wnere n g) an d \ n e ) are the eigenstates of 

the free Hamiltonian (corresponding to g — 0) and a n — arctan ( - ] . In this 



V S 

notation | n) are field number states and \g) and | e) are the ground and excited states of 
the qubit. We require that the qubit only acquires phases during its evolution. 
Therefore, to ensure that there are no transitions between the U(9) \n±) states, we 



make use of the adiabatic approximation 44 



|(m|H|n)| 



t < < 1. In the Jaynes- 
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Cummings model transitions are only possible between the eigenstates with same n, 

\{n + \H\n-)\ gt 

therefore the adiabatic condition yields 1 ^ — < < 1 . Notice that with 

' £„+-£„_ 2 

the values of g and f that we consider in the text the adiabatic condition holds. For 
simplicity, we assume that the parameter 6 undergoes a cycle, therefore the 
interaction time is t = 2tiI5. The dynamical phase y^ lt acquired by the state \n±) 
under this condition is given by 



TBn = 



-4g 2 n 



(9) 



The geometric phase y£ n defined by 31 : iy^ n — < U*doU> n + d6, yields 

Jo 

y£ n — 2tz^h — cos 2 ^ — j j and y^ n — In^n — sin 2 ■ We are interested in the case 

where the field is initially in a thermal state p F — p n \tif}(nf\, where 



Pn = {e-y{\-e-%F-. 



k B T 



In the mixed case, the geometric phase is given by 



Yq — Arg(^^~^ p„e( iJc; ' i ) j 45 ' 4& , while in the dynamical case, the phase is obtained 

through a simpler expression y^ — pn)'D>r We find, using a series expansion on 
g/S, that the geometric phase is given by Equation (2) and the dynamical phase by 
Equation (3). 

Quantum fisher information. The quantum Cramer-Rao bound 47 states that the 
error 3x in estimating a parameter x with M measurements on a state p x that depends 
on x, is bounded by: 



3x> - 



(10) 



F(p x ) being the Quantum Fisher Information of the state. This inequality assumes 
that input states are separable, however, it has been shown that entangled states 

improve the 1 j \fhA shot-noise scaling reaching the Heisenberg limit given by a II M 

scaling where M is the number of entangled particles. F{p x ) can be written as 48 : 
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l(m|g,p»| 2 

Pm+P„ 



(ID 



where p nun are the matrix elements of p x in an orthogonal basis in which p x is 
diagonal: p x — p n \ n) {n\ . The terms in which both p m and p n are zero in the sum 
in Eq. (11) must be excluded. 



1 . Pethick, C. J. & Smith, H. Bose Einstein Condensation in dilute gases {Cambridge 
University Press, Cambridge, 2004). 

2. Bruschi, D. E., Friis, N., Fuentes, I. & Weinfurtner, S. On the robustness of 
entanglement in analogue gravity systems. New J. Phys. 15, 113016 (2013). 

3. Wilson, C. M. et al. Observation of the dynamical Casimir effect in a 
superconducting circuit. Nature 479, 376 (2011). 

4. Jaskula, J.-C. et al. Acoustic Analog to the Dynamical Casimir Effect in a Bose- 
Einstein Condensate. Phys. Rev. Lett. 109, 220401 (2012). 

5. Ketterle, W. & Zwierlein, M. W. Making, probing and understanding Bose- 
Einstein Condensates, {in Ultracold Femi Gases, Proceedings of the International 
School of Physics Enrico Fermi, IOS Press Amsterdam, 2006). 

6. Szczepkowski, J. et al. Analysis and calibration of absorptive images of Bose- 
Einstein condensate at nonzero temperatures. Rev. Sci. Instrum. 80, 053103 

(2009) . 

7. Gerbier, F. et al. Experimental study of the thermodynamics of an interacting 
trapped Bose-Einstein condensed gas. Phys. Rev. A 70, 013607 (2004). 

8. El-Badry, A. M., Hassan, A. S. & Soliman, S. S. M. Thermodynamic properties of a 
condensed 39 K Bose gas in a harmonic trap. PhysicaB: Condensed Matter 410, 63 
(2013). 

9. Leanhardt, A. E. et al. Cooling Bose-Einstein Condensates Below 500 Picokelvin. 
Science 301, 1513 (2003). 

10. Meppelink, R., Rozendaal, R. A., Roller, R. A., Vogels, R. A. & van der Straten, P. 
Thermodynamics of Bose-Einstein-condensed clouds using phase-contrast 
imaging. Phys. Rev. A 81, 053632 (2010). 

11. Stenger, J. et al. Optically Confined Bose-Einstein Condensates. /. Low Temp. 
Phys. 113, 167 (1998), ISSN 0022-2291. 

12. Gati, R., Hemmerling, B., Foiling, J., Albiez, J. & Oberthaler, M.-K. Noise 
Thermometry with Two Weakly Coupled Bose-Einstein Condensates. Phys. Rev. 
Lett. 96, 130404 (2006). 

13. van Zoest, T. etal. Bose-Einstein Condensation in Microgravity. Science 328, 1540 

(2010) . 

14. Modugno, G. et al. Bose-Einstein Condensation in Microgravity. Science 294, 
1320 (2001). 

15. Modugno, G., Modugno, M., Riboli, F., Roati, G. & Inguscio, M. Two Atomic 
Species Superfluid. Phys. Rev. Lett. 89, 190404 (2002). 

16. Giinter, K., Stoferle, T., Moritz, T., Kohl, M. & Esslinger, T. Bose-Fermi Mixtures 
in a Three-Dimensional Optical Lattice. Phys. Rev. Lett. 96, 180402 (2006). 

1 7. Simoni, A. et al. Near-threshold model for ultracold KRb dimers from interisotope 
Feshbach spectroscopy. Phys. Rev. A 77, 052705 (2008). 

18. Tung, S.-K. et al. Ultracold mixtures of atomic 6Li and 133Cs with tunable 
interactions. Phys. Rev. A 87, 010702(R) (2013). 

19. Chin, C. et al. Precision Feshbach spectroscopy of ultracold Cs2. Phys. Rev. A 70, 
032701 (2004). 

20. Recati, A., Fedichev, P. O., Zwerger, W., von Delft, J. & Zoller, P. Atomic Quantum 
Dots Coupled to a Reservoir of a Superfluid Bose-Einstein Condensate. Phys. Rev. 
Lett. 94, 040404 (2005). 

21. Diener, R. B., Wu, B., Raizen, M. G. & Niu, Q. Quantum Tweezer for Atoms. Phys. 
Rev. Lett. 89, 070401 (2002). 

22. Bohi, P. et al. Coherent manipulation of Bose-Einstein condensates with state- 
dependent microwave potentials on an atom chip. Nature Phys. 5, 592 (2009). 

23. Schlosser, N., Reymond, G., Protsenko, I. & Grangier, P. Sub-poissonian loading 
of single atoms in a microscopic dipole trap. Nature 411, 1024 (2001). 

24. Chuu, C.-S. et al. Direct Observation of Sub-Poissonian Number Statistics in a 
Degenerate Bose Gas. Phys. Rev. Lett. 95, 260403 (2005). 

25. Kollath, C, Kohl, M. & Giamarchi, T. Scanning tunneling microscopy for 
ultracold atoms. Phys. Rev. A 76, 063602 (2007). 

26. Spethmann, N. et al. Dynamics of Single Neutral Impurity Atoms Immersed in an 
Ultracold Gas. Phys. Rev. Lett. 109, 235301 (2012). 

27. Steffen, A. et al. Digital atom interferometer with single particle control on a 
discretized space-time geometry. Proc. Natl. Acad. Sci. 109, 9770 (2012). 

28. Esteve, J., Gross, C, Weller, A., Giovanazzi, S. & Oberthaler, M. K. Squeezing and 
entanglement in a Bose-Einstein condensate. Nature 455, 1216 (2008). 



29. Ng, H. T. & Bose, S. Single -atom -aided probe of the decoherence of a Bose- 
Einstein condensate. Phys. Rev. A 78, 023610 (2008). 

30. Bruderer, M. & Jaksch, D. Probing BEC phase fluctuations with atomic quantum 
dots. New }. Phys. 8, 87 (2006). 

31. Fuentes-Guridi, I., Carollo, A., Bose, S. & Vedral, V. Vacuum Induced Spin-1/2 
Berry's Phase. Phys. Rev. Lett. 89, 220404 (2002). 

32. Fuentes-Guridi, I., Bose, S. & Vedral, V. Proposal for Measurement of Harmonic 
Oscillator Berry Phase in Ion Traps. Phys. Rev. Lett. 85, 5018 (2000). 

33. Martin -Martinez, E., Dragan, A., Mann, R. B. & Fuentes, I. Berry phase quantum 
thermometer. New }. Phys. 69, 053036 (2013). 

34. Bertlmann, R. A., Durstberger, K., Hasegawa, Y. & Hiesmayr, B. C. Berry phase in 
entangled systems: A proposed experiment with single neutrons. Phys. Rev. A 69, 
032112 (2004). 

35. Orth, P. P., Stanic, I. & Le Hur, K. Dissipative quantum Ising model in a cold-atom 
spin-boson mixture. Phys. Rev. A 77, 051601(R) (2008). 

36. Cirone, M. A., Chiara, G. D., Palma, G. M. & Recati, A. Collective decoherence of 
cold atoms coupled to a Bose-Einstein condensate. New /. Phys. 11, 103055 (2009). 

37. Hansel, W., Hommelhoff, P., Hansen, T. W. & Reichel, J. Bose-Einstein 
condensation on a microelectronic chip. Nature 413, 498 (2001). 

38. Meyrath, T. P., Schreck, F., Hanssen, J. L., Chu, C.-S. & Raizen, M. G. Bose- 
Einstein condensate in a box. Phys. Rev. A 71, 041604 (2005). 

39. Sinatra, A., Castin, Y. & Witkowska, E. Nondiffusive phase spreading of a Bose- 
Einstein condensate at finite temperature. Phys. Rev. A 75, 033616 (2007). 

40. Widera, A. et al. Entanglement Interferometry for Precision Measurement of 
Atomic Scattering Properties. Phys. Rev. Lett. 92, 160406 (2004). 

41. Posazhennikova, A., Birmuske, R., Bruderer, M. & Belzig, W. Entanglement 
generation in a system of two atomic quantum dots coupled to a pool of 
interacting bosons. Phys. Rev. A 88, 042302 (2013). 

42. Retzker, A., Cirac, J. I., Plenio, M. B. & Reznik, B. Methods for Detecting 
Acceleration Radiation in a Bose-Einstein Condensate. Phys. Rev. Lett. 101, 
110402 (2008). 

43. Martin-Martinez, E., Fuentes, I. & Mann, R. B. Using Berry's Phase to Detect the 
Unruh Effect at Lower Accelerations. Phys. Rev. Lett. 107, 131301 (2011). 

44. Ballentine, L. E. Quantum Mechanics: a modern development (World Scientific 
Publishing, Singapore, 1998). 

45. Sjoqvist, E. et al. Geometric Phases for Mixed States in Interferometry. Phys. Rev. 
Lett. 85, 2845 (2000). 

46. Singh, K., Tong, D. M., Basu, K., Chen, J. L. & Du, J. F. Geometric phases for 
nondegenerate and degenerate mixed states. Phys. Rev. A 67, 032106 (2003). 

47. Giovannetti, V., Lloyd, S. & Maccone, L. Advances in quantum metrology. Nature 
Photon. 5, 222 (2011). 

48. Paris, M. G. A. Quantum estimation for quantum technology. Int. J. Quant. Inf. 
07, 125 (2009). 



Acknowledgments 

We are indebted to Tomi Johnson for his very helpful insights. The authors would like also 
to thank Gerardo Adesso, Mehdi Ahmadi, Michael Berry, Per Delsing, Goran Johansson, 
Tim Ralph, Enrique Solano and Daniel Oi for useful discussions and the financial support of 
EPSRC Bridging the gaps. C.S. and I.F. acknowledge support from EPSRC (CAF Grant No. 
EP/G00496X/2 to I.F.) A.W. acknowledges funding from EPSRC grant No. EP/H027777/1. 

Author contributions 

I.F. conceived the main idea and directed the project, C.S. and AW. mainly executed the 
project, L.H. assisted in the experimental implementation. All the authors contributed to 
the discussion and manuscript preparation. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Sabin, C, White, A., Hackermuller, L. & Fuentes, I. Impurities as a 
quantum thermometer for a Bose-Einstein condensate. Set. Rep. 4, 6436; DOL10.1038/ 
srep06436 (2014). 



| This work is licensed under a Creative Commons Attribution 4.0 International 
License. The images or other third party material in this article are included in the 
article's Creative Commons license, unless indicated otherwise in the credit line; if 
the material is not included under the Creative Commons license, users will need 
to obtain permission from the license holder in order to reproduce the material. To 
view a copy of this license, visit http://creativecommons.Org/licenses/by/4.0/ 



SCIENTIFIC REPORTS | 4:6436 | DOI: 10.1038/srep06436 



6 



